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ABSTRACT 

We report on the analysis of all 65 pointed Rossi X-ray Timing Explorer observations of the recently 
discovered soft X-ray transient MAXI J1659-152 (initially referred to as GRB 100925A). The source 
was studied in terms of its evolution through the hardness-intensity diagram (HID) as well as its 
X-ray variability properties. MAXI J1659-152 traced out an anti-clockwise loop in the HID, which 
is commonly seen in transient low-mass X-ray binaries. The variability properties of the source, in 
particular the detection of type-B and type-C low-frequency quasi-periodic oscillations, and the way 
they evolve along the HID track, indicate that MAXI J1659-152 is a black hole candidate. The 
spectral and variability properties of MAXI J1659-152 imply that the source was observed in the 
hard and soft intermediate states during the RXTE observations, with several transitions between 
these two states. 

Subject headings: X-rays: binaries — X-rays: individual (MAXI J1659-152) 
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1. INTRODUCTION 

Black hole X-ray binaries have been studied since the 
early 1970s. It is now well established that the X-ray 
spectral and variability properties of these sources are 
strongly correlated, which is most clearly seen in the 
transient black hole X-ray binaries (BHTs) . While there 
is a great variety in the observed outburst behavior of 
BHTs (even for single sources), their outbursts typically 
proceed along 'q'-shaped tr acks in hardness-intensit y di- 
agrams (HIDs; s ee, e.g., iHoman fc Bellonil 120051 and 
iDunn et al.l [20Tol ). which are traced out in an anti- 
clockwise manner; we note that similar t racks are also 
traced out by neutron star transients (e.g. iTudose et al.l 
[200llLinaresl[20 09'). The various branches of these tracks 
correspond to distinct spectral states. Not all sources 
show all possible states and the time they spend in each 
state can differ significantly. 

There are various conventions for describing the spec- 
tral (aiid_jwiabilitYO s tates of BHTs; see, e.g., reviews 
bv IHoman fc Bellonil (|2005[ ). iRemillard fc McClintocH 
mm. In this p aper we follow the convention used 
in Bclloni (2010), which is based on the work by 
[Belloni et al. (2005) on GX 339-4. This source is often 
used as a template for the outburst evolution of BHTs, 
owing to its well defined q-shaped HID tracks, and the 
fact that it shows behavior that is common to many other 
systems. 

When GX 339-4 goes into outburst, its intensity is 
low and spectrum is hard (low- hard state or LHS). As 
the intensity increases, the spectrum remains hard, until 
the source makes a transition to the intermediate state 
(IMS), where the hard color starts to decrease at a rather 
constant intensity. The IMS can be divided into a hard 
and a soft IMS (HIMS and SIMS, respectively) depend- 
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ing on the spectral and variability characteristics ob- 
served. The transition from the LHS is always first to the 
HIMS. GX 339-4 often shows several transitions between 
HIMS and SIMS until the hardness decreases even fur- 
ther and the source reaches the so called high-soft state 
(HSS), where subsequently the hardness remains approx- 
imately constant as the intensity eventually decreases. 
At some point during this decrease hardness increases 
again and the source transits from the HSS via the IMS 
back to the LHS, and returns to quiescence. 

A BHT exhibits various types of quasi-periodic oscilla- 
tions (QPOs) and broad-band noise components, whose 
properties are strongly correlated with spectral state 
(see, e.g.. Klein- Wolt fc van der KlisH2008l ). The power 
spectrum of the LHS is characterized by strong broad- 
band variability (0.01-100 Hz fractional rms amplitude 
up to 50%). During t he HIMS so-called type-C QPOs 
are observed (see Wiina nds et al . 1999; Remillard et al.l 
120021 : iCasella et al.l 2005. for QPO type definitions), of- 
ten with strong harmonic content and accompanied by 
strong broadband variability (fractional rms up to 30%). 
In the SIMS various types of variability are observed: 
power spectra with type-B QPOs (also with strong har- 
monic content) or (weaker) type-A QPOs, but also power 
spectra with weak peaked noise and/or QPO features 
that have been poorly characterized. All SIMS power 
spectra have in common that their broadband variabil- 
ity is (considerably) weaker than in the LHS and HIMS. 
Type-C QPOs are typically observed between 0.1-10 Hz, 
whereas type A and B QPOs are generally confined to 
the 4-8 Hz range. In the HSS variability reaches its mini- 
mum strength (a few percent rms) ; sometirnes ve ry weak 
QPOs above 10 Hz are seen ()Homan et al.ll200"ll) . 

In this paper we present a study of Rossi X- 
ray Timing Explorer (RXTE) observations of MAXI 
J1659-152 (henceforth J1659). This source was dis- 
covered with th e Swif t Burst Alert Telescope (BAT; 
iBarthelmv et al.l 120051 ) on September 25, 2010 and 
was initially thought t o be a gam ma-ray burst (GRB 
100925 A, see Mangano et al.l [20101 ). Later, it was sug- 
gested by IKanni (j2010D to be a new Galactic X-ray 
transient due to its persistent X-ray emission; this was 
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confirmed by iNegoro et al.l ()2010D with MAXI obser- 
yations. J1659 has also been det ected in the radio 
(|van der Horst et alJl2010aD. opt ical fJelinek et alJ l2010D 
and sub-mm bands ()de Ugarte P ostigo ct al. 2010). The 
aim of this letter is to discuss the results of the aperiodic 
timing analysis and color evolution from RXTE observa- 
tions, based on which we conclude th at J1659 is a black 
hole candidate ()Kalamkar et al.ll2010[ ). 

2. OBSERVATIONS AND DATA ANALYSIS 

We analyzed all 65 RXT E Proportional Counter Array 
fPCA I Jahoda et all [200 6') observations of J1659 taken 
between September 28, 2010 (MJD 55467) and Novem- 
ber 8, 2010 (MJD 55508), after which the source was 
not observed due to solar viewing constraints. The in- 
tensity and colors were obtained using the Standard 
2 data from all the PCUs active during the observa- 
tions. The intensity is defined as the count rates in 
the 2.0-20.0 keV band; hard color as the count rate ra- 
tio of the 16.0-20 keV and 2.0-6.0 keV bands. Count 
rates are background corrected and normalized using the 
Crab nebula observations. This method of normaliza- 
tion (jKuulkers et al.l [19941 ) is based on the assumption 
that the Crab nebula is constant in intensity and hard 
color. It is used to correct for diff erences between PCUs 
(see, e.g. Ivan Straaten et al.ll20Q5l for details). 

The timing analysis was done with the Event mode 
data in the 2-60 keV range. Power spectra were gener- 
ated for each observation by averaging the fast Fourier 
transform power spectra of continuous 128-s intervals, 
with a Nyquist frequency of 4096 Hz and a lowest fre- 
quency of 7.8 mHz. Periods of dipping activity in th e 
X-ray light curves, as reported bv lKuulkers et al.l (|2010l ). 
were not excluded from our timing analysis. The Pois- 
son nois e spectrum is estimated with the analytical func- 
tion of iZhang et al.l (|1995D and subtracted from this 
average power spectrum. The resulting power spec- 
trum is expressed i n source fractional rms normalization 
(Ivan der KlisllT989l ). using the average background rate 
dm-ing the observation. Each power spectrum was fit 
with a multi-Lorentzian (1-5 components) function and 
(when necessary) a power law. 

We performed a spectral analysis for only a few se- 
lected observations, including those cases where HID lo- 
cation and timing properties were inconclusive with re- 
spect to the spectral state of the source. Spectra were 
extracted from the Standard 2 mode data of PCU2, using 
HEASOFT v6.10. The spectra were background sub- 
tracted, corrected for dead time, and a systematic er- 
ror of 0. 6% was appli ed. Fits were made with XSPEC 
V12.6.0 ()Arnaudlll99"6l ) between 3.0-40.0 keV, with phe- 
nomenological models consisting of a simple accretion 
disk component (diskbb), a (cut-off) power-law, and a 
Gaussian line (fixed at 6.4 keV). Reported fiuxes are in 
the 2-20 keV band and are corrected for interstellar ab- 
sorption, which was fixed to 1.7 x 10^^ a toms cm~^, the 
avera ge value in the direction of J1659 (jKalberla et al.l 
l2005l) . 
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3. RESULTS 

color evolution 



Light curves ana 

The top panel of Figure [T] shows the 1.5-20 keV light 
curve from MAXI. It includes the early rise and late 
decay, which were not covered by our RXTE data set 
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Fig. 1. — Top to bottom: From MAX/ observations - (A) in- 
tensity (cnts cm~^ s~^, 1.5-20.0 keV); from iJXTi? observations - 
(B) intensity (2-20 keV), (C) hard color (16.0-20.0/2.0-6.0 count 
rate ratio), (D) average 0.01-100 Hz fractional rms amplitude, and 
(E) QPO frequency. Each point represents one observation, except 
for the MJD 55488.02 and 55490.11 observations, for which we use 
two points each. Different types of power spectra are indicated 
using different symbols as follows: (i) open circles: type-C QPOs 
which indicate the HIMS; (ii) filled triangles: type-B QPOs; (ill) 
filled circles and filled squares are observations with weak variabil- 
ity (iv) stars are the observations with unclassified power spectra. 
The latter three categories comprise the SIMS. 

(see below). Panels B-E of Figure [D show the RXTE 
2-20 keV light curve, hardness curve, time evolution of 
fractional rms amplitude (0.01-100 Hz), and QPO fre- 
quency evolution, respectively. The outburst showed 
a relatively rapid rise; when it was first observed with 
RXTE, three days after discovery, it already had a flux 
of 6.6xl0~® ergs/cm^/s. The rise was followed by an 
extended plateau during which some irregular intensity 
variations (time scale of a few days) were observed (here- 
after referred to as 'flares'). During two observations in 
this flaring phase (MJD 55488.02 and 55490.11) we ob- 
served several rapid ('^20 s), almost step-function-like 
changes in the count rate from the source. The changes 
were on the order of 30%, lasted ^100-1500 s, and 
were not accompanied by changes in the spectral hard- 
ness. Similar phenomena have been observed in, e.g. , 
GX 339-4 and XT E J1859-h226 ()Mivamoto et al.l 119911: 
iCasella et al.l I2004D and have been referred to as 'flip- 
flops'. For both observations the high and low count 
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Fig. 2. — With the RXTE data, we plot A) Hardness-intensity 
diagram (HID) and, B) 0.01-100 Hz fractional rms amplitude vs. 
hard color. Each point represents one observation, except for the 
MJD 55488.02 and 55490.11 observations, for which we use two 
points each (see text). Different power spectral types are indicated 
(symbols as in Figure [T]!. The locations of three of the power 
spectra shown in Figure[3l fPDSl. PD2, and PDS5) are indicated 
with arrows. With the MAX/ data, C) HID, the intensity and hard 
color measured in 1.5-20.0 keV band and in 10.0-20.0/2.0-4.0 keV 
band, respectively. 



rate levels were separated and treated as different obser- 
vations for the remainder of our analysis (and are plotted 
separately in Figures [TH2|) . The flux reached its maxi- 
mum during a flare on October 8 (MJD 55477) at a value 
of lxlO~^ ergs/cm^/s . A sudden radio flux decre ase was 
reported for that day (jvan der Horst eralll2010b[ ). From 
the intensity plateau the source evolved towards a slow 
decay, after which it presumably returned to quiescence. 
The flux during the last i^XTi? observation was 1.4x 10~^ 
ergs/cm^/s. 

Figure [2j\ shows the HID. The source started in the 
upper-right corner and traversed its track in an anti- 
clockwise manner, ending in the lower-right corner. Hys- 
teresis is observed in the spectral evolution, with the 
hard-to-soft spectral evolution occurring at a count rate 
^3 times higher than the soft-to-hard spectral evolu- 
tion. While the spectral evolution during the early and 
late phases of the outburst was largely monotonic, the 
source exhibited frequent back-and-forth motion in spec- 
tral hardness during the flaring phase of the outburst. 

Since RXTE did not observe J1659 during the 
first three days of its outburst, we used MAXI 
(jMatsuoka et al.|[2009() data to infer the spectral evolu- 
tion during the early rise. We plot the MAXI HID in 
Figure [2lC. Overall, the path traced out is similar to that 
seen in the RXTE HID. The first part of the MAXI HID 
(i.e. the right vertical branch and subsequent left turn) 
was not covered by RXTE. A comparison with HIDs of 
other BHTs suggests that during the first three days of 
the outburst, J1659 evolved through parts of the LHS 



a nd HIMS. 

iShaposhnikov fc Yamaokal ()201Cl[ ) reported that during 
the softest part of the outburst J1659 reached the HSS. 
We performed a spectral fit of the observation with the 
lowest spectra hardness (MJD 55489.26) to verify this. 
We obtain an excellent fit (reduced x^=0.96, for 62 d.o.f.) 
with a 0.8 keV disk black-body and a power-law compo- 
nent with index 2.2 (plus a weak Gaussian at 6.4 keV). 
The disk componen t contributed ~60% to the 2-20 keV 
flux. iBellonil ()2010[ ) does not provide a spectral defini- 
tion of the HSS, but following the definition of the soft 
or thermal dominant) state bv lRemillard fc McClintoc^ 
2006[) . which uses 75% as a lower limit, we conclude that 
the source did not reach the soft state during this out- 
burst. This behavior is similar to, e.g., the 2000 outburst 
of XTE J1550-564 (Miller ct al. 2001). 



3.2. Variability properties 

Figure [2j3 shows the relation between hard color and 
broad-band rms. Such a diagram can be helpful in identi- 
fying different type s of power spectra ((Belloni et al.|[2005l : 
iFender et al.l 120091) . At high hard colors, hardness and 
rms are correlated, but at low hard colors (<0.5) there 
is considerable scatter. 

Based on the overall strength of the broad-band vari- 
ability, the spectral hardness, and the shape of the power 
spectra (noise/QPOs) we divided our RXTE observa- 
tions into five groups (see symbols in Figures [T][2]) . For 
the classification of the PDS and the QPOs, we ni a de use 
of the cl assificatio n schemes of iRemillard et al.l (|2002[ ) 
and Ca seUa et al.l ([2005). Representative power spectra 
are shown in Figure [3] and are referred to as PDS 1-5. 

3.2.1. Group 1: type-C QPOs 

The Group 1 PDS have strong 0.01-100 Hz variability 
(~23%), a strong QPO (varying between 1.6-8.5 Hz), 
at times with a harmonic and a sub-harmonic. They 
were observed during the spectrally hard parts of the 
outburst (hard color > 0.45; see open circles in Figures 
[T] and . The frequency is strongly anti-correlated with 
spectral hardness (Fig.[TJi;). The QPO frequencies and Q 
values, and the fractional rms amplitude in the 0.01-100 
Hz frequency range are typical of PDS exhibiting type-C 
QPOs, as observed in the HIMS (no phase lags or energy 
dependence of the QPOs were measured). Hence, we 
identify these QPOs as type-C. 

3.2.2. Group 2: type-B QPOs 

The Group 2 PDS have 0.01-100 Hz variabflity < 10% 
rms and power-law noise as opposed to the strong vari- 
ability and peaked noise of the Group 1 PDS. These PDS 
are observed in a fairly narrow range of spectral hardness, 
have QPOs with frequencies between 2.6 and 4.1 Hz, at 
times accompanied by a harmonic and a sub-harmonic, 
and cluster in a small patch of points in the hardness-rms 
diagram. In Figure|4l we plot the 0.01-100 Hz fractional 
rms amplitude vs. frequency of these and the PDS with 
type-C QPOs, as well as the PDS with the unclassi- 
fied QPOs (see below). Wha t we observe i n Figu re IH is 
very similar to the results of iCasella et al.l (|2004D . indi- 
cating that these are type-B QPOs, which classifies these 
observations as SIMS. Fractional rms amplitude anti- 
correlates with frequency for type-C QPO, while type-B 
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Fig. 3.— Five representative power spectra, labeled PDS 1-5. PDS 1 (ObsID: 95108-01-02-00) and PDS 2 (ObsID: 95118-01-01-01) are 
typical of type-C and type-B QPOs, respectively. PDS 3 is an average power spectra of the five observations classified as 'weakl'. PDS 4 
is the average power spectra of the nine observations classified as 'weak2'. PDS 5 (ObsID: 95118-01-12-00) is the power spectrum of an 
observation which had an unclassified QPO. The HID locations of PDS 1, 2 and 5 are indicated in Figure[2] 



QPOs do not follow such an anti-correlation. Note that 
after the detection of the first two type-B QPOs (at 3.9 
Hz and 4.1 Hz), the spectrum temporarily hardened dur- 
ing the next two observations, and type-C QPOs (at 6.0 
and 6.9 Hz) were observed, (see Fig.[T]) followed again by 
type-B QPOs (3.9 and 3.5 Hz). This indicat es rapid 
SIMS/HIMS transitions ()Wiinands et al.l Il999f ) over a 
period of ^^6 days. 

3.2.3. Groups 3 & 4' weak power 

In the soft part of the outburst most power spectra 
do not show clear indications for QPOs. Such power 
spectra were observed during two time intervals. During 
the first time interval (MJD 55486-55590) the source in- 
tensity moved up and down on a daily basis — this is 
also the period during which the flip-flops were observed 
(see i jS.ip . Inspection of the individual power spectra 
revealed subtle differences between the low count rate 
and high count rate observations. PDS3 is the averaged 
power spectrum of the five high count rate level obser- 
vations (weakl in Figures [T] and [2]) ; it is consistent with 
a single power-law with no distinguishable features. The 



low-count rate observations show indications for an ad- 
ditional broad bump around 7-8 Hz. Such a feature is 
also seen for the observations in the second interval with 
weak power (MJD 55496.5-55498.5). PDS4 is the aver- 
aged power spectrum of these nine observations (weak2 
in Figures [T] and [2]) . The PDS4-like power spectra oc- 
curred at slightly lower hardness than the PDS3 ones. 
Power spectra similar to PDS3 an d PDS4 were observed 
in GX 339-4 (jBelloni et al.l 120051 ) and were grouped to- 
gether with power spectra that show type-A (not ob- 
served here, see below) and type-B QPO to define the 
SIMS. Therefore, also based on this definition J1659 did 
not reach the HSS. 

3.2.4. Group 5: unclassified QPO 

For three power spectra that showed somewhat broad 
QPOs (Q~1.6-3.5) we had difhculty assigning a classi- 
fication — they are marked as unclassified in Figures 
[iHH The observations were among the shortest (800 s) 
and the power spectra were of lower quality. However, 
the locations of these observations in Figure |4] (and also 
Fig. [2j3) suggest that these QPOs may be type-B QPOs 
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Fig. 4. — 0.01-100 Hz fractional rms amplitude vs. frequency 
of the strongest QPO in each power spectrum. The type-B and 
typc-C QPOs form different branches, with the unclassified QPOs 
falling close to the type-B ones. 

with a lower frequency (^1.9 Hz) and also lower Q- value, 
which would classify these observations as SIMS. 

We did not see indications for type-A QPOs in our 
observations. These are typically accompanied by simi- 
larly shaped noise as type-B QPOs, but show no strong 
harmonic content (as seen in PDS2). Also in the rms- 
frequcncy plot they fall in a d ifferent location than type- 
B QPOs (.Casella eran[200a . which argues against the 
QPOs in the type PDS5-like power spectra being of type- 
A. 

3.2.5. Search for High-frequency QPOs 

High-frequency QPOs in the the 50-450 Hz 
range have been report e d for a handful of BHTs 
(jMcClintock &: RemillardI I2OO60 . We inspected aU 
individual power spectra in the total 2-60 keV band, as 
well as in the 5.7-60 keV band, but no indications for 
such QPOs were found. We also inspected the averaged 
power spectra of the different groups (with groups 2 and 
5 added together). Again, no significant QPOs were 
found, with a 2.1a measurement of a peak at ^160 Hz 
(Q~4, ^4% rms) in the 2-60 keV power spectrum of 
PDS 4 constituting our most significant detection. 

4. DISCUSSION & CONCLUSIONS 

We analyzed all the RXTE observations of 
MAXI J1659-152 taken during its 2010 outburst. 
Spectral states were identified using the HID in combi- 
nation with the aperiodic variability properties. In the 
HID J1659 traced out an anti-clockwise loop, similar to 
those seen in black hole and neutron-star transients (see 
iJIJ. The types of power spectral features we identified 



(in particular type-B and type-C QPOs) and the way 
they evolve along the HID track, indicates that J1659 is 
a black hole candidate. The spectral decomposition of 
the softest spectrum (i.e. diskbb + po wer-law, withou t 
the need for second thermal component: ILin et al.l l2007l. 
the radio loudness in the early part of the outburst 
(jParagi et al.l I2010D . the fiip-flops seen in two obser- 
vations, and the sudden (~day) drop in the rms (i.e. 
PDS3) at lo w-hardness, a r e all characteristic of th i s type 
of systems (|Dunn et al.l 120101: IMivamoto et al.l 119911: 
IFender et al.ll2009D . strengthening our identification of 
J1659 as a black hole candidate. 

Based on spectral and variability properties we find 
that J1659 moved from a HIMS to a SIMS and back 
again (at lower intensity), without reaching the HSS 
(or the soft or thermal domi nant state as defined by 
'RemiUard iiMcClmtock'2006). Our i^XTO observations 
started after the source had left the LHS (as indicated by 
AL4X/ observations), and ended before it had returned 
to that state. 

The observations in the softest part of the SIMS had 
weak variability and revealed two types of power spectra, 
one consistent with a power-law and one with an addi- 
tional bump around 7-8 Hz. We note that similar groups 
of observations with w eak variability were o bserved in the 
SIMS of GX 339-4 bv lBelloni et "all ([2001 . IFender et al.l 
(|2009f) found that radio flares and subsequent quenching 
of the radio fiux often occur in a time interval of a few 
days before and after the time of sudden drops in the 
rapid X-ray variability (identified as a distinct zone in 
their rms-hardness diagrams). In J1659 a similar drop 
occurred on MJD 55481, a few da ys after the quench- 
ing of t he radio flux as reported by Ivan der Horst et al.l 
( j2010b[ ), consistent with other BHTs. 

Optical observations constraining the mass of the com- 
pact object are required to confirm the black hole na- 
ture of J1659. From dipping episodes in the X-ray light 
curves, the or bital period period ha, s been proposed to b e 
2.4-2.5 hour (jKuulkers et al.|[20lol: iBelloni et afllMol) . 
This makes the source very interesting as, if confirmed, 
this would make it the black hole binary with the shortest 
known orbital period. 
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